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ABSTRACT 

We present high-resolution, 77-band, imaging observations, collected with Subaru/HiCIAO, of the 
scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub- 
mm imagery suggested the existence of the dust-depleted cavity at r < 46 AU, our observations reveal 
the presence of scattered light components as close as O'/ 2 (~ 28 AU) from the star. Moreover, we 
have discovered two small-scale spiral structures lying within O'/ 5 (~ 70 AU). We present models for 
the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h ~ 0.1, 
which is consistent with previous sub-mm observations. This model can potentially give estimates of 
the temperature and rotation profiles of the disk based on dynamical processes, independently from 
sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable 
on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely 
identify the origin of these spirals, planets embedded in the disk may be capable of exciting the 
observed morphology. Assuming that this is the case, we can make predictions on the locations and, 
possibly, the masses of the unseen planets. Such planets may be detected by future multi- wavelengths 
observations. 

Subject headings: cireumstellar matter — instrumentation: high angular resolution — polarization 
— protoplanetary disks — stars: individual (SAO 206462, HD 135344B) — waves 
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Based on data, collected at the Subaru Telescope, which is operated by the National Astronomical Obser\atory of Japan. 
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1. INTRODUCTION 

Dynamical processes in protoplanetary disks such as turbulence or disk-planet interaction are important in under- 
standing physical condition and evolution of disks, and planet formation processes. High resolution, direct imaging 
observations of circumstellar/protoplanetarv disks can reveal non-axisymmetric structures, providing insight into such 
dynamical processes (e.g., Hashimoto et al. 2011). 

Recent observations have identified a class of protoplanetary disks harboring tens of AU-scale holes/gaps at their 
centers: the so-called transitional disks. One well-studied system of that class is the rapidlv rotating Bering F star, 
SAO 206462 (HD 135344B, F4Ve, d = 142± 27 pc, M = 1.7 tg;?/U 0 , Muller et al. 2011). The observations of CO line 
profiles (Dent et al. 2005; Pontoppidan et al. 2008; Lyo et al. 2011) and stellar rotation (Muller et al. 2011) consistently 
indicate an almost face-on geometry ( i 11°). The gap in the disk was predicted from the infrared (IR) spectral 
energy distribution (SED, Brown et al. 2007), and was subsequently imaged in sub-mm dust continuum at ~ O'/ 5 x0"25 
resolution (Brown et al. 2009). Andrews et al. (2011) estimate the gap radius to be ~ 46 AU and the surface density 
within the gap to be 10 -52 times smaller than that extrapolated from the outer disk. The gas in Keplerian motion 
surrounding the gap region is also imaged by CO lines (Lvo et al. 2011). The CO rovibrational line observations 
(Pontoppidan et al. 2008) and [OI] spectral line observations (van der Plas et al. 2008) indicate the presence of a 
gas disk in the vicinity (several AU-scale) of the star. SED modeling (Grady et al. 2009) and NIR interferometry 
(Fedele et al. 2008) indicate the presence of an inner dust belt, which is temporally variable (Sitko et al. 2012) and not 
coplanar with the outer disk (Benisty 2011, private communication). New imaging with high spatial resolution and 
sensitivity is required to understand the inner structures of the disk. The outer portions of gaps can now be resolved 
using 3-10 m ground-based telescopes at near infrared (NIR) wavelengths (e.g., Thalmann et al. 2010, for LkCa 15). 

In this Letter, we present 77-band polarized intensity (PI) observations of the disk of SAO 206462 down to r ~ 
O'/ 2 (~ 28 AU) scale at Of' 06 (~ 8.4 AU) resolution. Interior to the sub-mm resolved gap, we find spiral structures, 
indicative of dynamical processes. We use the spiral densitv wave theory to interpret the structure, and estimate disk’s 
physical parameters. 


2. OBSERVATIONS & DATA REDUCTION 

2.1. HiCIAO Observations 

SAC* 206462 was observed in the 77-band (1.6 pm) using the high-contrast imaging instrument HiCIAO (Tamura et al. 
2006; Hodapp et al. 2008; Suzuki et al. 2010) on the Subaru Telescope on 2011 May 20 UT as part of Strategic Ex- 
plorations of Exoplanets and Disks with Subaru (SEEDS, Tamura 2009). The adaptive optics system (AO 188; 
Hayano et al. 2004; Minowa et al. 2010) provided a stable steliar point spread function (PSF, FWHM = 0'.'06). We 
used a combined angular differential imaging (ADI) and polarization differentia! imaging (PDI) mode with a field of 
view of 10 ’ x 20 ' and a pixel scale of 9.5 mas pixel -1 . A O'.'3-diameter circular occulting mask was used to suppress 
the bright stellar halo. The half-wave plates were placed to four angular positions from 0°, 45°, 22.5°, and 67.5° in 
sequence with one 30 sec exposure per wave plate position. The total integration time of the PI image was 780 sec 
after removing low quality images with large FWHMs by careful inspections of the stellar PSF. 

2.2. PDI Data Reduction 

The raw images were corrected using IRAF 1 for dark current and flat-field following the standard reduction scheme. 
We applied a distortion correction using globular cluster M5 data taken within a few days, using IRAF packages 
GEOMAP and GEOTRAN. Stokes (Q, U) parameters and the PI image were obtained in the standard approach 
(e.g., Hinkley et al. 2009) as follows. By subtracting two images of extraordinary- and ordinary-rays at each wave 
plate position, we obtained +Q, — Q, +U, and —U images, from which 2 Q and 2U i mages were made bv another 
subtraction to eliminate remaining aberration. PI was then given by PI = \/Q 2 + U 2 . Instrumental polarization of 
HiCIAO at the Nasmyth instrument was corrected by following Joos et al. (2008). From frame-by-frame deviations, 
the typical error of surface brightness (SB) was estimated to be ~ 5% at r ~ O'/S when averaged over 5x5 pixels (~ 
PSF scale). Comparing different data reduction methods (frame selections and instrumental polarization estimates), 
we expect that the systematic uncertainty of the SB of PI to be ~ 10%. 

2.3. Contemporaneous Photometry 

Since SAO 206462 shows variability in NIR wavelengths (Sitko et al. 2012), it is important to take photometry 
simultaneously with disk observations. 77-band photometry was obtained just before and after the disk imaging 
without the coronagraphic spot with the adaptive optics, by sixteen 1.5 sec exposures at four spatially dithered 
positions. An ND10 filter (9.8 ± 0.1% transmission) was used to avoid saturation. Using the MKO filter set, the 
77-band (A e s = 1.615 pm, FWHM = 0.29 pm, Tokunaga et al. 2002) magnitude was 6.96 ± 0.07 mag. 

Broad-band VRIJHK photometry was obtained on 2011 May 23-26, starting within 48 hours of the HiCIAO ob- 
servation, using the Rapid Eye Mount (REM) Telescope at La Silla, Chile (Covino et al. 2004). The REM 77-band 
filter has A e fr = 1 .65 pm, FWHM = 0.35 pm: broader and displaced to longer wavelengths than the MKO filter. The 
observed data were reduced differentially using SAO 206463 (A0V). The IR excess due to the inner disk (Figure 1) was 

1 IRAF is distributed by the National Optical Astronomy Obser-atory, which is operated by the Association of Universities for Research 

in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Fig. I.— SED for SAO 206462 obtained by REM, adapted, in part, from Sitko et al. (2012). The REM observations consist of VRIJHK 
photometry. Also shown are spectra obtained with the SpeX spectrograph. See Section 2.3 for the data reduction techniques. 


average for the range observed in 2009-2011 (Sitko et al. 2012). No significant variation was observed during the 2011 
May observations, except for the small long-term fading trend (5m = 0.08 ± 0.02 mag) over the observation period. 

Figure 1 also displays spectra obtained with the SpeX spectrograph (Ravner et al. 2009) on NASA’s Infrared Tele- 
scope Facilitv (IRTF). The observations were obtained in the cross-dispersed (XD) echelie mode between 0.8 and 5.1 /mi 
using a 078 slit ( R ~ 900) and calibrated using HD 129685 (A0V) with SpeXtool (Vacca et al. 2003; Cushing et al. 
2004) . The absolute flux calibration, to correct for light loss at the spectrograph slits, was accomplished in two ways: 
using photometry and wide-slit spectroscopy (see Sitko et al. 2012). The March SpeX data were normalized using the 
REM photometry, obtained in the days immediately after the SpeX observations, and when the star was photometri- 
cally stable. In Julv, SAO 206462 and the calibration star were observed with the SpeX prism using a 3'/0 slit, which, 
under good seeing and transparency conditions, produces absolute flux to ~ 5% accuracy. The REM photometry at 
H- band in May is ~ 0.25 mag brighter than the Subaru data. The Subaru photometry shows a low value even when 
considering the filter difference, suggesting that the outer disk is illuminated efficiently. 

3. RESULTS 
3.1. Spiral Structure 

The SAO 206462 disk can be traced in PI from O'/ 2 to ~ l'/O (28 — 140 AU), similar to the range imaged with 
HST/NICMOS (Grady et al. 2009), but with a factor of ~ 4 greater angular resolution. The total PI is 9.87 mJy ± 
0.06% at O'/ 2 < r < l'/0, which is 0.6% of the stellar intensity. The total PI at 0742 ' r < 170 is 3.94 mJy ± 0.1% 
while the total intensity by HST/NICMOS F110W is 9.7 mJy (Grady et al. 2009). The average SB of PI at r = 0746 
is ~ 6 mJy/asec 2 , whereas the total intensity by HST/NICMOS F160W is 30 mjy/asec 2 (Grady et al. 2009). Given 
the NICMOS data uncertainties, the polarization fraction is ~ 20 — 40%, assuming no PSF halo in the HiCIAO data 
and no variable self-shadowing/illumination in the disk. Our measured polarization fraction is similar to that of HD 
100546 (14%^g^, Quanz et al. 2011) and AB Aur (~ 25% — 45%, Perrin et al. 2009). 

Figure 2 shows the PI image. The region interior to 074 is not a void and we do not see clear structural evidence 
of the cavity wall in Andrews et al. (2011) model (R cav = 46 AU ~ 0733). We see spiral arcs SI (east) and S2 (south- 
west). The PI at the location of the spirals is ~ 30% larger than that extrapolated from the smooth outer profile 
(bottom of Figure 2). The brightest portions of the spirals roughly coincide with the bright thermal emission peaks at 
12 pm (Marinas et al. 2011) and lie inside the ring noted by Doucet et al. (2006). It is also noted that we see a dip in 
PI in the north-west, probably due to the depolarization in the minor axis direction (see below), and that we do not 
see large-scale, localized shadow that might be cast by the inner dust belt if highly inclined relative to the outer disk. 
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Fig. 2. — Top: PI image of SAO 206462 in the north-up configuration with log-stretch color scales. The filled orange circles at the center 
indicate the mask size (r = 0"15). The circles around have r = O'.' 2, exterior to which the features are considered to be real. The right 
panel is central region's closeup. Different color scales are used to enhance the spirals labeled as SI” and “S2” . The “Dip” may be due to 
depolarization. Bottom: PI profile along the red line in the top right panel. The arrow indicates the location of S2. The position errors 
are not shown for visibility. 


3.2. Azimuthal and Radial Profiles 

Here, we summarize the overall disk structure exterior to the spirals. Figure 3 (top panels) shows the azimuthal SB 
profiles at r = C//5 and r = OUT. At r £ Of! 5, SB has maxima around position angle (PA) ~ 50° — 60° and 230° — 240°. 
Since the polarization is maximized at ~ 90° scattering (e.g., Graham et al. 2007), it is implied that the disk major 
axis is at PA ~ 50° — 60°, comparable to estimates by CO observations: PA = 56° ± 2° by Pontoppidan et ai. (2008) 
and 64° ± 2° by Lyo et al. (2011). We adopt PA = 55° for the major axis and i = 11° for inclination (see Section 1). 
Our spiral model fitting results (the next section) are little affected even if we assume a face-on geometry. 

From CO observations, it is known that the south-west side is receding (e.g., Lyo et al. 2011). Therefore, either the 
north-west or south-east side is the near side. We do not see an obvious forward scattering excess in the NIR image. 
However, since the spirals are typically trailing, it is inferred that the south-east is the near side. 

Figure 3 also shows the radial PI profiles along the major axis, which is roughlv consistent with r -3 in the outer 
part, indicating a flat (not highlv flared) disk (Whitney and Hartmann 1992). The radial slopes vary as PA from ~ —2 
to ~ -4.5 (fitting at Of.'Q < r < 1"0) or from ~ —2.5 to ~ — 5 (fitting at O'/ 3 < r < Of! 9), with shallower slopes typically 
appearing in the minor axis directions. However, r~ 3 is representative on average. This slope is observed in several 
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Fig. 3. — Top: Azimuthal SB profile at r = (X'5 (ieft) and r = C//7 (right) with PA measured from north to east. Bottom: Radial profile 
along the major axis. Position errors indicate the FWHM of the PSF. SB errors are estimated using frame-by-frame deviations. 

other HAeBe disks (e.g., Fukagawa et al. 2010, for total intensity data), although HD 97048 disk exhibits a shallower 
slope (Quanz et al. 2012). 


4. SPIRAL STRUCTURE MODELING 

Among several features in the image, the most interesting one is the rion-axisymmetric spirals. In order to understand 
them, we propose a model based on the spiral densitv wave theory (e.g., Lin and Shu 1964; Goldreich and Tremaine 
1978, 1979; Ogilvie and Lubov 2002), assuming that NIR emission traces the disk surface density structure. With 
such a model, the spiral structures can be used to infer the disk temperature, independently of, for example, CO line 
observations. 

The shape of the spiral density wave is determined by the location of the launching point (corotation radius r c ) and 
disk’s thermal and rotation profiles. When the disk rotation angular frequency is fl(r) oc r~ a and the sound speed is 
c(r) cc r~P, the shape of the wave far from r c is given by 

9(r)= e 0 _5£fcn) * 

fllQ 


7nV 4 

»f 1 

1 ( r Y 

a \ ( 1 i Al 

\r c J 

\l + /3 

1 -a + 6 \r c J 

j \l + 0 1 — q + /V 


in the polar coordinate (r, 9), where h c = c(r c )/r c fl(r c ) denotes the disk aspect ratio at r = r c and 8q gives the phase. 
Equation (1) well approximates the shape of the density wave given by the WKB theory (Rafikov 2002; Muto et al. 
2011). When the spiral is excited by a planet in a circular orbit, its location is ~ (r c , # 0 ). Equation (1) has five 
parameters, (r c ,9o,h c .a,/3). 

Two non-axisymmetric features, SI and S2, (Figure 4) are identified as follows. First, local maxima in the radial SB 
profiles normalized by r 2 (to take into account the dilution of the stellar flux) are traced at every 1° step with data 
at 170° < PA < 360° (SI) and at 50° < PA < 190° (S2). The points near the minor axis (PA ~ 325°) are excluded 
because the structure there may be affected by depolarization. The points at PA > 200° (Si) and PA > 114° (S2) 
may be a part of axisvmmetric rings since they are found to have constant radii. After removing these points, we have 
27 (SI) and 56 (S2) points as representing samples of non-axisymmetric spirals, with the opening angle of ~ 15° for 
both 31 and S2. We estimate that the uncertainty of the location of the maxima is given by the FWHM of the PSF. 
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Fig. 4. Left: Red and green points indicate the locations of the maxima of radial profiles for SI (top) and S2 (bottom). The green ones 
are used for the fitting by Equation (1) (Section 4). Dashed lines indicate circles with r = 0" 39 (top) and r — Cr/3 (bottom). Right: The 
spiral shape given by Equation (1) with the best-fit parameters of SI (top) and S2 (bottom). Crosses show the locations of (r c ,8 o). 


In order to fit the non-axisymmetric structures by Equation (1), we fix a and 8 at 1.5 (Kepler rotation) and 0.4, 
respectively, as in Lyo et al. (2011), while other parameters are varied as (O'/I < r c < O'/ 9, 0 < 9q < 2-rr, 0.05 < h c < 
0.25). Note that different values of /3 yield similar results. Since it is difficult to fit SI and S2 simultaneously, they are 
fitted independently. 

The ‘ best-fit ,! parameters are (r c ,6o,h c ) = (0"39, 204°, 0.08) for SI (reduced y 2 = 0.52) and (r c , 9q, h c ) = 
(0'.'9, 353°, 0.24) for S2 (reduced x 2 = 0.31). The spiral shapes with these parameters are shown in Figure 4. However, 
the parameter degeneracy is significant. Figure 5 shows the parameter space of (r c ,0 o ) with 63.8% confidence level 
for h c = 0.1 and h c = 0.2. Note that in Figure 5, the ‘’best-fit ,: of (r c , 0o) is outside the domain of confidence in some 
cases because h c is not the same as the best-fit. Despite the parameter degeneracy, the values of the aspect ratio 
which fit the shape of the spiral ( h c ~ 0.1) are consistent with those obtained from the sub-mm map of the disk (e.g., 
h = 0.096(r/100AU)° 15 ; Andrews et al. 2011). 

The spiral density wave theory predicts that the pattern speed deviates from the local Kepler speed; 


n 


pattern 



[deg/vr] 


( 2 ) 
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FIG. 5.- Parameter degeneracy of the fitting for SI (top) and S2 (bottom) with = 0.1 (left) and h c = 0.2 (right). The red points 
show a part of the data fitted. The green curves indicate the 63.8% confidence level for the locations of (r c ,0o). The magenta crosses are 
(r c , &o ! for the best-lit parameters. Note that the best-fit values of h c are different from those shown in the figure. 

is not necessarily equal to 0(r). When r,. = 0/5 (<•'' 70AU), the spiral will move ~ 10° in a decade, corresponding to 
movement of O'/I. Considering the PSF scale of our observations and the locations of the spirals, such deviations can be 
detectable over a couple of decades. Moreover, if the two spirals have distinct corotation radii, their relative locations 
change in time due to the pattern speed difference. Such measurements will confirm that the observed feature is really 
the density wave, providing indisputable evidence of dynamical activity. 

Note that it would be difficult to detect spirals in colder disks (smaller h c ), where spirals are more tightly wound, 
due to the blurring by the PSF. The lower detectable limit of h c is typically h c ~ 0.01 — 0.03 for our set of parameters. 
The combination of high angular resolution and warm temperatures allows the spiral structure in the SAO 206462 
disk to be resolved. Further spirals might be detectable in similarly warm disks. 

5. SUMMARY AND DISCUSSION 

In this Letter, we present a high-resolution image of the SAO 206462 transitional disk using Subaru/HiCIAO, with 
an inner working angle of 0/2. We discover non-axisymmetric spiral features, which can be explained bv the spiral 
density wave theory with a reasonable value of the disk aspect ratio ( h c ~ 0.1). The model is robust in a sense that 
it does not assume the origin of such structure explicitly. 

The detection of scattered light within the sub-mm cavity itself is interesting, since Andrews et al. (2011) predicts 
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that the sub-mm cavity is heavily depleted. Our data in tandem with the millimeter data may suggest that the 
depletion of grains at different sizes is not uniform. Dong et al. (2012) discuss such discrepancies between sub-mm 
dust continuum images and NIR scattered light images in terms of differing spatial distributions as a function of grain 
size from a general theoretical perspective. 

Our major assumption is that PI at if- band, tracing the scattered light at the disk surface, actually traces the surface 
density variations. This assumption is valid when the disk is in vertical, isothermal hydrostatic equilibrium without 
rapid radial surface density variations (e.g., Muto 2011). Structures near the midplane are, however, preferentially 
observed at longer wavelengths with high spatial resolution: Atacama Large Millimeter /Submillimeter Array (ALMA) 
can be an ideal instrument. 

Among several possible causes for the spiral structures (see also Hashimoto et al. 2011), one interesting idea is that 
planets excite them. In this case, the domain of possible locations of the planets is given by the green curves in Figure 
5. If the two spirals have distinct corotation radii, there may be two (unseen) planets embedded in the disk. The 
amplitude of the surface density perturbation scales with the planet mass as <fE/E ~ GAfpO/c 3 for non-gap-opening 
low-mass planets ( GM p Q/c 3 1, e.g., Tanaka et al. 2002). In our data, the amplitude of the spiral wave is typically 

~ 30% (Figure 2), implying M p ~ 0.5 Mj. 

The typical error of PI in our image is ~ 5%. If this is typical of HiCIAO, it is capable of detecting the indirect 
signatures of planets down to M p ~ 0.05Mj. ADI is promising in finding a point source; however, small field rotation 
(~ 13 3 ) due to the southerly declination (—39°) of SAO 206462 makes obtaining the total intensity difficult in our 
data. Future /^band observations may reveal thermal emission from a planet, if it exists, or its surrounding (accreting) 
gas. 
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